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Competitive exclusion, a key principle of ecology, can be generalized to understand many other 
complex systems. Individuals under surviving pressure tend to be different from others, and correla- 
tions among them change correspondingly to the updating of their states. We show with numerical 
simulation that these aptitudes can contribute to group formation or speciation in social fields. 
Moreover, they can lead to power-law topological correlations of complex networks. By coupling 
updating states of nodes with variation of connections in a network, structural properties with 
power-laws and functions like multifractality, spontaneous ranking and evolutionary branching of 
node states can emerge out simultaneously from the present self-organized model of coevolutionary 
process. 

PACS numbers: 89.75.Hc, 89.20.Hh 



Process of competitive exclusion[l[ occurs in some real 
systems-evolutionary branching of species in ecosystems, 
citations in scientific research and designation of con- 
sumer goods are examples among many others. It is 
actually a fundamental ingredient governing main prop- 
erty of dynamical behaviors of systems which are often 
described with complex networks^ nowadays. However, 
the contribution of competitive exclusion to the interac- 
tional structure of networks and to their functional fea- 
tures is not widely realized up till now. 

In modeling a system, individuals are represented as 
nodes and correlations among them are represented as 
edges of a graph. Scale- free property Q, characterized 
by power-law degree distribution, has attracted exten- 
sive attention since it reflects a general feature of di- 
verse systems such as the Internet, citation networks, 
protein-protein interaction, and so on[2j. In most pre- 
vious models, dynamics of networks and dynamics on 
the networks are separated. The interplay between the 
formation of topological structure and functions emerge 
from the network is usually neglected, which is reasonable 
when the structure is independent of the dynamical states 
of nodes, or when these two sides vary in rather differ- 
ent speeds. However, in many practical phenomena like 
academic and art creation, financial transactions, global 
climate fluctuation and synaptic plasticity of neuron net- 
work in the brain both the structure and functions 
emerge from the identical process, and time-dependent 
variations of both individual states and local connections 
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of nodes feed back with each other. Therefore, novel 
models with coevolution mechanisms Q underlying them 
appeared to fit for the necessary. Unfortunately, scarcely 
could one produce both scale-free structure and collec- 
tive dynamics of nodes simultaneously. On the other 
side, new nodes are often assumed to know the global 
information of whole the growing network, which is usu- 
ally impossible for huge-size systems. In this sense it is 
needed to set up models based on local interactions to 
see if structure and functions at system level will emerge 
from self-organized dynamics 

As it is well known, competitive exclusion plays key 
role in the formation of species. There is strong com- 
petition among species occupying the same or nearest 
loci, surviving pressure force them drift their traits away 
from the local average level, and gradually induces evo- 
lutionary branching of species. Sympatric speciationQ 
in an ecosystem is a recent focus of naturalists. It refers 
to the origin of two or more species from a single local 
population. Seceder model[l] based on a simple rule of 
local third-order collision succeeded in mimicking such a 
process and capturing its similarity to group formation in 
society. A network version Q of it has been reported, giv- 
ing rise to a possible mechanism of community structure 
and clustering in social networks. 

In this paper the principle of competitive exclusion 
is generalized outside the realm of ecology, the seceder 
model is modified to describe temporally updated states 
of nodes and corresponding variation of connections 
among them together. We show that generic natures of 
members in diverse systems, i.e., to be different from oth- 
ers under the pressure of competition, and coevolution 
between updating node states and varying connection 
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among nodes, can lead to simultaneous emergences of 
evolutionary branching of individual traits, spontaneous 
ranking and multifractality of node states and, power-law 
topological structure of correlations in a system. In this 
way we are able to understand scale-free phenomena and 
other characteristics in various fields with a novel com- 
mon mechanism. Such self-organized coevolution mod- 
els of scale-free network with both structural and func- 
tional properties integrated are still few to the best of 
our knowledge. 

We set up the present model through three iteration 
rules. (1) Network growth starts from a primitive com- 
plete graph with too nodes. Each node i on joining the 
network was assigned an initial state with a random real 
number w(i) uniformly distributed in (0, 1). At each time 
step, a new node i' is added to the preexisting network. 
It gives out to edges (to < m ) to old nodes arbitrarily. 
(2) At every step, each node i counts w(i) -the average 
of state values w(j)(j ^ i) over its nearest linked neigh- 
bors, from them it picks up the one whose w(J) makes the 
maximum distance from average w(i), i.e. J m ax{i) cor- 
responds to max\w(j) — w(i)\, then, a randomly selected 
node j among the nearest neighbors of i is chosen as the 
offspring of Jmax(/t), called J se d{i)- Different from origi- 
nal seceder model|8( , it is kept at its own site and, with its 
state variable updated as w{J se d(i)) = w(J max (i)) + 5, 
where random number S € (0, 1) is also uniformly dis- 
tributed and with positive numerical range for wider ap- 
plications. Obviously w(i) here can be accounted as a 
time-dependent non-decreasing fitness^. (3) For the 
new comer node i' at every step, together with its 'young' 
enough fellows (i.e. %' — i < AI, with AI a given inte- 
ger constant implicating aging effect fill] . hereafter we call 
them I altogether for convenience). Search seceders for 
all I's neighbors j . When w(J se d(j))/w(I) > h, where h is 
a given value of threshold, a new edge is added between 
node J se d{j) and /(Double links and self-loops are for- 
bidden). Meanwhile, an edge linking such node / and its 
neighbor j is removed if the condition w(j)/w(I) < h or 
w(I) / w(J) < h is satisfied. Finally, if any node i becomes 
isolated due to edge-cutting, directly link it to its seceder 
J se d(i)- The threshold description of correlation adop ted 
here is widely used in modeling complex systems 
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Actually the iteration rules of the model are abstracted 
from observation to real systems. In art creation and sci- 
entific research, people have generic tendency to create 
new works so that they behave differently from others. 
Sparkles from collision of opinions with large difference 
often result in creation. As well known, scholars are of- 
ten under the pressure of publication. Papers with the 
same or very similar viewpoint, method and results to 
existing ones have less chance to get published. Here we 
see the competition exclusion promotes prosperity of sci- 
entific research. Suppose a graduate student just start 
his academic career by joining the research on a certain 
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FIG. 1: (color online) Power-law degree distributions of co- 
evolving scale-free networks, (a) Degree distributions with 
h = 3.0, lines for various values of m collapse onto a single 
one with 7 = 2.39; Inset: In-degree distribution with (3 = 2.0 
for h = 3.0 and m = 1. (b) Threshold-dependent degree dis- 
tributions with h — 3.0, 4.0, 5.0 and 6.0 respectively. Inset: 



Corresponding in-degree distributions, 
and mo = 20 for all lines. 



N = 10 4 , AI = 10 



topic, usually he has to focus on some papers after exten- 
sive searching due to limited time, and often he extends 
his reading to references of them. Generally speaking, 
he needs to pay more attention to ones with sharp con- 
trasts against his knowledge background(w(i)), and un- 
derstand lately published literature (w(j sed(j))) to in- 
spire new ideas for his own paper. But in the reading he 
may be restrained within the ability of his understand- 
ing. Therefore, it is natural to predict a suitable range 
of threshold ratios within which papers with state values 
w(J sed(j)) would be cited (connected). And papers in 
selective reading based on one's local sight are likely to 
be cited, forming increased in-degree of those ones. On 
the opposite, papers(on the node state w(j) ) have small 
difference (too low ratio of w(j)/w(i)) with w(i) are less 
cited(the link between node i and j is trimmed). Any- 
way, a recently updated node state(w(jsed(i))) would be 
more attractive to a failure (an isolated node). Artists up- 
date themselves by continuous creation, therefor the co- 
occurrence network of musicians serves another example 
of competitive exclusion. We know that musicians with 
similar genre are competitors for performance. Managers 
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FIG. 2: (color online) (a)Asymptotic size-dependent power- 
law decay of Pearson coefficients: r(N) ~ N~ a for N > 10 3 , 
where a depends on thresholds h. (b) Size-dependent power- 
law decay of clustering coefficient, C(N) ~ N~ v . Averaged 
on 10 realizations of network configurations. 
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FIG. 3: (color online)For low valuesi.e./i = 2.0, dashed line, 
and h = 1.5, full line, respectively), degree distributions de- 
viate from power-law, but they are retrieved by introducing 
ten percent (circles for h = 2.0) and twenty percent (squares 
for h — 1.5) relaxatio respectively on carrying out rule (3). 
Other parameters are the same as those in Fig. la. 



usually do not intend to arrange opportunity for them to 
appear on the same stage since audience prefer perfor- 
mance with diversity. It is assumed that whoever created 
a playlist was using a certain criterion to group artists in 
them. One does not normally find concerts with a mix- 
ture of heavy rock, jazz and piano sonata, therefore a 
range of thresholds is used to balance the homogeneity 
and heterogeneity. As the results of coevolution, both ci- 
tation network [l3, EH EH an d musician network[l6j dis- 
play the topologyof scale-free structure although most 
foodwebs do not[I3|. Suppose a man faces to job crisis, 
he has to refresh himself to become non-trivial for going 
out of dilemma. And he may attempt to learn from, even 
coalesce to a succeeded person by recommendation of a 
common friend. But whether they can sustain a close re- 
lation, it depends on whether they are mutually needed 
and compensate in a proper measure (e.g. w(i)). In all 
these cases states of nodes keep varying with time and 
correlations among them change corresponding to such 
variations along an optimal gradient. 



Coevolution of node states and topological connec- 
tion yields most structural properties of complex net- 
works by self-organization. Numerical simulation reveals 
out power-law distribution of node degree: p(k) ~ fc~ 7 , 
which is illustrated in Fig. 1 a. Without ensemble aver- 
age on network configurations, it is shown that in the case 
of h = 3.0 the distribution keep invariant for all values 
of m, with the slope 7 = 2.39. In-degree is counted by a 
node to its accepted edges from younger ones. The dis- 
tribution also shows essential power-law as shown in the 
inset of Fig. la. The slope of the double-logarithmic line 
Pi(k) ~ fc _/3 is around [3 = 2.0, which is in accordance 
with numerical results of another model [l3| and empir- 
ical studies [l^. [TH . In Fig. lb wc show the variation of 
power exponents 7 depending on correlation thresholds 
h. They lay in the range of (2.0, 3.0), which fits well 
to real complex systems. And the inset of it displays 



that essential power-law behavior of in-degree distribu- 
tions also exist for different thresholds. The calculated 
Pearson coefficients r[l8| which describe degree-degree 
correlation of the network are shown in Fig. 2a. They are 
positive reflecting statistical feature of social networks. 
Moreover, they also behave asymptotic power-law decay 
in the size of the system, i.e. r(N) ~ N~ a , which is, 
to our knowledge, a specific feature and first predicted 
by the present model. It is expected to be verified by 
empirical data from real complex systems. Fig. 2b dis- 
plays size-dependent decay of clustering coefficients 0: 
C(N) ~ N~ v . The exponents arc (1.2 — 1.4) correspond- 
ing to thresholds in the range of [3.0,8.0] while for ran- 
dom graphs we have 77 = 1 for comparison. 

Simulations with low threshold values (e.g., h = 2.0 
and 1.5) reveal some different behaviors of the coevolu- 
tion, where iteration rules no longer lead to power-law 
degree distributions, (see full and dash lines in Fig. 3). 
However, when we allow a small portion (ten to twenty 
percent) of cut-off operations not to carry rule 3, scale- 
free properties can be retrieved promptly (sec Fig. 3). 
Moreover, degree-degree correlations restore assortativ- 
ity corresponding to it. This implies that randomness 
may play an essential role in the origin of scale-free be- 
haviors since there should be more or less relaxation on 
deterministic rules in complex systems (l9j|. 



Ranking behavior of node states also emerge spon- 
taneously from coevolution. Whole the range of node 
states is divided into 100 intervals in Fig. 4 to show that 
the values are distributed quite discontinuously. This 
is drastically different from uniform initial distribution 
and, is comparable to group formation in original seceder 
model(see Fig.l of ref.Q). Inherited from seceder model, 
two prominent traits(see Fig. 4) at both ends can be re- 
garded as the result of evolutionary branching Q with 
the tendency of elimination for mediate genotypes. Here, 
species in sympatry seem to likely drift their traits away 
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FIG. 4: Histogram for numbers of nodes N(L) with discrete 
ranks L of node state values w(i), h — 3.0 and m = 18, other 
parameters are the same as in Fig. la. 
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FIG. 5: Check the existence of state multifractality of node ac- 
cording to consecutive joining the network with box-counting 
method. Inset: Singularity spectrum /(a) of multifractal of 
node state w(i) for the network with h = 3.0, mo = 20 and 
m = 18. 



from local average level since the- strongest competition 
exists between similar genotypes [2(|. Anyway, to make a 
scrutiny into applicability of co-evolutionary mechanism 
to sympatric speciation would be valuable. Applied to 
citation networks, it means that the long term coevolu- 
tion gradually eliminate the publishing chance of a paper 
at middle level, instead, the population of quality tends 
to be divided and shift approaching both ends. Beyond 
seceder model]!, @, our numerical results also give sup- 
port to the assumption of the ranking model [2l[ of SFN 
with self-organization mechanism. It is noticeable, scale- 
free property as a result of coevolution can be obtained 
without the prerequisite of preferential attachment on the 
power-law function of prestige ranks of nodes. 



The updating process of node states induced by com- 
petitive exclusion coupled with topological variation re- 



sults in collective behavior of nodes, which reflects char- 
acteristics of functional aspects apart from structural 
ones of the network. Based on simulated data of node 
states w(i) which are put in order of the time sequence 
as node's participation in the network, we calculate func- 
tion V(q,d) = y), d)lnni(q, d) with standard box- 
counting technique [23 for different moment q versus 
x =ln d, where d represents scales of boxes, and fii is 
normalized measure of the summation over states in box 
I. Essential linearity can be seen for at least 4-5 center 
lines in Fig. 5 so that the verified singularity spectrum 
f(a) of multifractal is shown in its inset. Interestingly, 
the present work gives another example of long-range cor- 
related gradient-driven growth of a multifractal entity p3j 
with scale-free network as its inherent skeleton. The mul- 
tifractality of node states is found to emerge accompanied 
with scale- free property of the structure and vanishes cor- 
respondingly. We have verified the correspondence be- 
tween two properties in the range of too £ [15,50], and 
AI E [5,15]. Therefore, the present model suggests a 
common mechanism of scale-free structure of social sys- 
tems together with their multifractality and assortativity 
as well. 

Simultaneous emergences of macroscopic properties on 
both structural and functional sides also enable us to 
understand functions in coordination with the Internet, 
world wide spatial distribution of population^ with all 
kinds of transport and communication networks connect- 
ing resident sites being complex networks among which 
some ones are coevolution SFNs, middle latitude cli- 
mate network [12). citation network [l3l [Til , [poll, number 
distribution of species in ecological networks [24| , mu- 
sician networks jig, and diversity maintenance method 
for evolutionary optimization algorithms [l5l [25| . on a 
novel platform of coevolution with alterable details. Ac- 
tually, it relics on the mechanism with another type 
of preferential attachment of node state correlation but 
does not explicitly depend on node degree 0, [H, HH, 
which distinguishes itself from previous ones. Starting 
from but outreaching seceder model, we can account 
generic natures of individuals-to update states to self- 
adapt the competitive exclusion, and correlations among 
them change correspondingly-as driving force in self- 
organization of some evolutionary complex systems char- 
acterized by power-law distributions of various topologi- 
cal quantities and specific functions. 
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